Identification of transverse spin currents in noncollinear magnetic structures 
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We show that the transverse components of spin current in a ferromagnet is linked to an oflt diag- 
onal spin component of the transmission matrix at interfaces;it has little to do with the mismatch of 
band structures between dissimilar metals. When we take account of this component, not considered 
in prior analyses, we find spin torque comes from a region of at lease 3 nm around an interface. 
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Attention has recently been focused on current induced 
switching of magnetic layers; this idea was proposed in 
199611 and has been verified in the past three years. 
It occurs when one drives an electric current across a 
multilayered structure in which the magnetic layers are 
noncollinear. The central idea behind the switching is 
that when the spin polarized current that develops in a 
fixed ferromagnetic layer of a multilayered structure im- 
pinges on a second noncollinear free magnetic layer the 
component of the spin current transverse to the magne- 
tization is absorbed. As one assumes the conservation of 
spin angular momentum this creates the torque on the 
background magnetization to switch it. One of the unre- 
solved issues is the length scale over which the transverse 
component of the spin current is absorbed and what ma- 
terial parameters controls it. Here, we point out, that 
the detail of the transverse spin near the interface does 
matter in determining the spin transfer torque. The rea- 
son is that transport is inherently non-local so that the 
incident spin polarization at an interface does depend on 
the detail of the transverse spin accumulation even if this 
length scale is small, i.e., one must include transverse pin 
accumulation. Mathematically, the boundary condition 
solely from the vanishing transverse spin density in fer- 
romagnets is not sufficient to determine the spin torque. 
If we write the boundary condition to the distribution 
function on the two sides of an interface, we have found 
that the truncation of the transmission matrix to a diago- 
nal matrix in spin space misses key relations between the 
spinor distribution functions on the two sides of an inter- 
face, and thus the transport across the entire structure 
is unable to be self-consistently determined. 

We recently found a possible mechanism by which one 
can inject off diagonal spin distributions into each sheet 
of the spin split Fermi surface of a magnetic layer, i.e., 
in the presence a current a spin flip potential can ex- 
ist at interfaces in noncollinear magnetic multilayers. 



This additional scattering excites the transverse compo- 
nents of spin currents in magnetic layers that are found 
by adopting the Boltzmann equation's definition of spin 
current and then we find there is no discontinuity in the 
spin current at an interface between two dissimilar layers. 
This transverse current is different from that evaluated 
from equilibrium states; 0,0] they complement one an- 
other but the existence of the latter does not negate the 
former. We evaluate the relevant parameters over which 
this component is absorbed by using ab-initio band struc- 
ture calculations and find the characteristic length scale 
in a typical 3d transition-metal ferromagnet (Co) is 3 nm. 
This is an order of magnitude greater than that found in 
prior analysis which evaluated a quantum mechanical ex- 
pectation value of the current over the Fermi surface. 0,0] 
As the layers being switched are typically 2-3 nm this dif- 
ferences is important as it alters one's estimates of the 
critical current for switching. 

The spin current density at a point r can be defined 
asli 



Jcpir) = — (*t(^)v,f^(r) 



(1) 



where the expectation value of the operator is taken over 
a nonequilibrium distribution function and a, (3 are spin 
indices. This "quantum mechanical" definition of spin 
current Eq.QJ has been evaluated for a selected set of 
equilibrium states so as to make a spin polarized current 
enter a ferromagnet at an angle to the magnetization 
rather than self-consistently determing the nonequilib- 
rium distribution. 0,01 Based on this definition the trans- 
verse component of the spin current has been shown to be 
carried by the correlation between the states of opposite 
spin on the Fermi surface Um = kp] , km = kp^ . |j] This 
identification has lead one to conclude that it is the band 
mismatch between the electronic structures of adjacent 
layers that causes the polarization of the spin currents to 
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change abruptly in a narrow region about an interfaces 
so that the component of the spin current transverse to 
the magnetization of a layer is absorbed within about a 
i nm of the interface0,0|, i.e., the spin split band struc- 
ture of the magnetic layers does not allow a spin current 
transverse to the magnetization to propagate past the 
interfacial region. 0, Q We recently showed how in the 
presence of a current one can inject a spin distribution 
transverse to the magnetization across an interface with 
a well defined momentum; |^ here we determine how it 
propagates in a homogenous layer. 

The equation of motion for the electron transport den- 
sity matrix {o^ka^k'i3^, i.e., the quantum or semiclassical 
Boltzmann equation, in a translationally invariant sys- 
tem and for linear response is an equation for a given 
momentum state in which states at other momenta enter 
through the scattering terms. The Wigner transform 
of these matrix elements. [loj 
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is independent of position r for a system in equilibrium, 
e.g., the Fermi distribution function. While the pertur- 
bation of the electric field in and of itself does not break 
the translational invariance of a layer the ensuing cur- 
rent crossing an interface does; therefore when a current 
is driven across a system as a multilayer the out of equi- 
librium part of fafsik, r) develops a position dependence. 
In the Boltzmann description the transverse components 
of the spin current are 
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where in linear response v{k) is the Fermi velocity of the 
electron. In principle both definitions Eqs.lQJ and Q 
give the same current provided sufficient care is taken in 
their evaluation. 

The equation of motion for the density matrices or dis- 
tribution functions is the Boltzmann equation; here we 
use the semiclassical version of this equation. In each 
layer the energy and density vary slowly on the length 
scale of the Fermi wavelength so that we can limit our- 
selves to the first term in the gradient expansion of the 
equation of motion. When we limit ourselves to linear 
response in the external electric field we find the equa- 
tions of motion in steady state for the elements of the 
spinor density matrix for each momentum state on the 
Fermi surface kp are.[inl| 

v;d,fp-eEvpS{e~SF)^-^^^^^-i^^^^, (4) 



and 



(5) 



where we have used a simplified index p to denote the 
momentum kp of a state on the n*'' sheet of the Fermi 
surface(we suppress this index), p' are states of opposite 
spin to p, the average (/p) represents elastic scattering to 
all states on the Fermi surface, Vp is the component of the 
Fermi velocity along the electric field E, J is the magnetic 
part of the energy (see below) and we have limited the 
current induced variations of the distribution function to 
those along the growth direction of a multilayered struc- 
ture X (also the field direction). The diagonal elements 
fp = faa{kp, x) represents the occupancy of the state kp; 
in equilibrium it is given by the Fermi function so that 
only the spin state a that crosses the Fermi level is oc- 
cupied while the other is zero and we do not consider 
it further. The off diagonal elements fp ^ fuikp.x), 
which we call a current induced spin coherence, |9| repre- 
sent coherences between the state kp on the Fermi sur- 
face and the states with opposite spin; see Eq.lQl); these 
coherences occur when we drive a spin current across a 
normal metal/ferromagnet (N/F) interface. Q The scat- 
tering terms include those for states of the same spin on 
the Fermi surface as well as those between sheets 
of opposite spin r^^^; as we include the lat- 

ter only to have well defined boundary conditions on our 
distribution functions. 

The solutions for the longitudinal components are well 
know.fl^ From Eqs.Q and © we see that the elec- 
tric field only creates out of equilibrium longitudinal 
components of the distribution functions; in a homoge- 
neous magnetic layer there is no coupling to the trans- 
verse components. However when the spin current from 
one layer is injected into another noncoUinear magnetic 
layer the transverse components can be excited. As 
such they have to be determined by a self consistent 
calculation of the distribution functions in the differ- 
ent layers of a multilayered structure. When we ne- 
glect collisions (the rhs of Eq.©) the transverse solu- 
tions are f^{x) ~ exp±i{Jp/hVp)x; when we average 
this over the Fermi surface we find an interference be- 
tween individual p or k states so that the transverse 
components of the spin currents in the magnetic layers, 
jxi^) ~ / Vx{k)f'^{k,x)dk (Eq.lPJ), can be fit to a form 
approximating an exponential decay ^ exp — x/Afr-fl^ 
In the ballistic regime Xtr = dj = hvp/J where the 
bar denotes an average over states on the n*'' sheet 
of the Fermi surface under consideration: [l^ while for 
diffusive systems where we consider the collision terms 
Atr — Xj = -v/ d TX„, fp /3TT. Therefore when compo- 
nents of the spin current are injected into a magnetic 
layer that are transverse to its magnetization, we find 
they propagate a distance Xtr before decaying; as this 
distance is an order of magnitude greater than the Fermi 
wavelength one can describe the transverse spin currents 
in the semiclassical Boltzmann approach. As Xtr is com- 
parable to the thickness of the magnetic layers undergo- 
ing switching one cannot assume the transverse compo- 
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nents of spin currents are entirely absorbed in such thin 
layers. 

To evaluate dj we need the Fermi velocity w^(k), and 
J(k) for each point k on the n'^ sheet of the Fermi sur- 
face under consideration; the latter is identified as fol- 
lows. The energy of a spin polarized band in a uniform 
ferromagnetic metal can be written as a spinor 

en{k)=en{k)l + ]^Jn{k)M-a (6) 

where n is the band index, M is a unit vector in the di- 
rection of the magnetization, and a represents the Pauli 
spin matrices. As such J„(fc) is the magnetic part of the 
band energy. For a band crossing the Fermi surface it 
is found by identifying the band of opposite spin that 
is split from the common energy £„(fc) by J„(fc). This 
requires one to do two band structure calculations, one 
nonmagnetic in which each band is spin degenerate and 
another for the fully spin polarized state in which the 
two spin split states originating from the common de- 
generate band can be easily identified along directions of 
high symmetry. With this identification we find 

J„(fc) = (e„T(fc)-£«i(fc)), (7) 

where £na{k) are the energies of the spin split band and 
n is the index for the band that crosses the Fermi surface 
whose magnetic component we are evaluating. Therefore 
Jn(k) can be found from looking at the band structure 
of a ferromagnetic metal; it is the difference in energy at 
constant k between two split bands. While these two 
bands mostly lie on top of (next to) one another, there 
are regions of the Brillouin zone where one has to do some 
further analysis in order to identify the bands which had 
a common origin in the nonmagnetic phase. 

From the band structure for fee Co we find dj for Ai 
majority states along the T — X (100) direction is 7.4 nm 
but it decreases rapidly off this axis, along T — K (110) dj 
is 2.93 nm, along F — i (111) is l.lnm, and the average 
dj for the majority band over the Fermi surface is about 
3 nm.|l7j Of the three minority spin bands crossing the 
Fermi surface the one that is the partner of the majority, 
#3, has a dj of 0.8 nm while the other two have even 
shorter dj. Central to the semiclassical description of 
the transverse spin currents is the condition that they 
can propagate over distances large compared to the Fermi 
wavelength; we find that this condition is well satisfied 
for the 3d transition metal ferromagnets for the majority 
band. 

To identify the source of the discontinuity in a spin 
current at interfaces we have solved the Boltzmann equa- 
tions Eqs.Q) and lO for the current across two magnetic 
layers; while there is an intervening nonmagnetic layer we 
adopt the usual assumption that its thickness is much less 
than the spin diffusion length so the spin current is con- 
stant across it. For the electronic structure we take a free 



electron-like model in which the degenerate spin bands 
are split by the internal fields; the effect of the exchange 
splitting of the bands is to produce two separate chan- 
nels of conduction on the Fermi surface, one majority 
one minority. ^3 To keep track of current induced coher- 
ences the conduction in each channel is described by a 
2x2 density matrix so that one has to match two 2x2 
density matrices at interfaces to obtain a description of 
the spin current across a noncollinear multilayered struc- 
ture, i.e., each spin channel is represented by a spinor 
distribution function. To perform a self consistent calcu- 
lation of the distribution functions across a multilayered 
structure we have to match them at the interfaces by 
using transmission and reflection coefficients. In addi- 
tion to the spin diagonal coefficients Tmm^ss'{s' — s) 
previously considered we add the off diagonal coefficient 
with s' ^ s that transmits spin coherence across an N/F 
interface; Q therefore we must consider in our match- 
ing not only states on the Fermi surface but also those 
with which electrons form current induced coherences in 
steady state, i.e., in general we have to match 4x4 rather 
than 2x2 distribution functions across interfaces. 

When one limits oneself to spin diagonal transmission 
coefficients at an interface there is no transmission into 
off diagonal or transverse elements of the density matrix, 
i.e., the fn, are zero (not excited) in the magnetic layers; 
as they have been identified as the transverse components 
of the spin current one arrives at the conclusion that the 
transverse component of the spin current has been ab- 
sorbed at the interface. 0, Q However, when we include 
the off diagonal transmission coefficient the discontinuity 
of the transverse component of the spin current at inter- 
faces is removed. In other words without Tmm^^ss' with 
s' ^ s we find the spin current in a ferromagnetic layer 
is parallel to the background magnetization and there is 
an abrupt rotation of the spin current when the layers 
are noncollinear; with it the spin current is continuous. 
While transverse spin currents do not exist when one uses 
the equilibrium distribution that was selected to evalu- 
ate Eq.c3), y, 13 they do exist when we use Eq.© to 

self-consistently evaluate the transverse component of the 
spin current Eq.JSJ in steady state. |0| When we solve the 
Boltzmann for the spin current across two noncollinear 
magnetic layers we find its polarization starts to devi- 
ate from the background magnetization a distance Xtr 
on either side of the interface and is continuous across 
it.^2Q] While the spin accumulation, and therefore the 
resistance, depends on the reflection and transmission co- 
efficients due to the specular scattering at the interface, 
the continuity of the spin current is insensitive to them 
as long as we include T„„„^^ss'with s' 7^ s ; therefore 
we conclude that their discontinuity is not an indigenous 
property of noncollinear magnetic layers. We find it is 
caused when one suppresses the transverse components 
of the spin currents in the magnetic layers either by us- 
ing only diagonal transmission coefficients ^ or by using 
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only states at the Fermi surface to evaluate the spin cur- 
rent Eq.dl.li,!! 

In summary, when the steady state spin current is self- 
consistently evaluated using Eq.© we keep track of the 
transverse spin accumulation around the interfaces irre- 
spective of its length scale. We find there is no absorp- 
tion at the interface, it is continuous across it, and the 
transverse component decays within 3 nm, at least for 
Co, irrespective of the reflection and transmission coef- 
ficients; albeit these coefficients affect the magnitude of 
the spin currents and resulting torques. Therefore the 
discontinuity in spin currents at interfaces is not due to 
band mismatch; rather we find it arises from omitting the 
off diagonal components of the transmission coefficients 
and the ensuing transverse spin accumulation. The dif- 
ference between our result and those which predict spin 
torques created within ^ nm of an interfaces is impor- 
tant as the thickness of layers undergoing magnetization 
reversal is in the range of 2 — 3 nm. When we introduce 
diffuse scattering in the layers we reduce the decay length 
for the transverse component of the spin current as we 
found in the spin diffusion model of current induced mag- 
netic switching. [2 L] but we do not create discontinuities 
in the current. 
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